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I. INTRODUCTION 


The use of radar for the remote quantitative measurement of rain 
rate provides a powerful tool for those concerned with the distri- 
bution and occurrence of rainfall in the atmosphere. A knowledge 
of the spatial and temporal distributions of rainfall is essential 
in such diverse disciplines as meteorology, hydrology, and milli- 
meter wavelength communications. Radar has been extensively used 
for applications such as these since the development of meteorological 
radar techniques in the 1940's. Nevertheless, the use of this approach 
has been handicapped by a lack of adequate precision in the resulting 
quantitative rain rate data: for example, Harrold and Nicholass [1] 
report average uncertainties in such radar measurements ranging from 
30 to 100». These uncertainties are largely attributable to the 
difficulties encountered in performing an adequate calibration of 
the radar measurement system. 

Traditionally two distinct techniques have been used for radar 
calibration. The first method utilizes a metallic sphere or other 
target having a known radar cross section to establish an absolute 
calibration. Unfortunately, this method is based upon the use of a 
point target and does not compensate for the effects of beam-filled 
targets such as those usually encountered in rain rate measurements. 
Furthermore, it is often quite difficult to place such a standard 
calibration target in the far iield of the radar antenna and at a 
height sufficient to eliminate ground reflections. The second common 
method is to calibrate the radar observed rain rate against the rain 
rate measured on the ground below a given radar resolution cell. 

This approach obviously ignores differences resulting from both 
vertical and horizontal wind components since the two measurements 
to be compared are made at two different spatial locations. This dif- 
ficulty is compounded by the fact that the radar measurement must 
be elevated in order to avoid ground reflections. 

This report addresses the problem of radar calibration from 
several different directions. It is shown in the following that 
the integrated, range-corrected radar backscatter saturates in a 
manner analogous to that found in radiometric temperature measure- 
ments. Furthermore, it is shown that this saturation leads to an 
absolute calibration method for radars operating at attenuating wave- 
lengths. The method for calibration of meteorological radars in terms 
of independent radiometric temperature measurements is also developed. 
This technique offers the advantage that the calibration is based 
upon an independent measurement made throughout the same spatial 
region as that viewed by the radar. 
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II. THE RADAR RANGE EQUATION AND ITS SOLUTION 

The development below follows closely that originally presented 
by Hitschfeld and Borden [2]. Here it is assumed that a narrow 
beam, short pulse radar is viewing a region of rainfall which fills 
each radar resolution cell homogeneously. Then, if the power pattern 
of the radar antenna is assumed to be Gaussian, the following radar 
range equation results [3]: 

(1) Pr = ^ Z e RD [watts] 

r • 


where P r is the time average received power backscattered from a 
resolution cell located at a range of r meters from the radar. The 
radar calibration constant, Q, is 


m n - * 5 k 2 f 2 Gt Pt 

(2) Q * — [watts/meter] 

210 • tn 2 • c 

where 

iki 2 - 

n = the complex refractive index of liquid water 

(3) f - radar frequency [H z ] 

G « maximum antenna gain 

t a radar pulse length [sec.] 

Pt = peak transmitted power [watts] 

c = velocity of propagation [m/sec]. 

The equivalent reflectivity, Z, has been defined in terms of the 
volumetric radar cross section, [m2/m3]: 


n -1 


n2+? 


(4) 


Z = 


c 4 n 


,5r4 


f 4 |k| 


[m®/m3]. 


Ard is the one-way attenuation between the radar and the resolution cell 
being observed expressed in terms of specific attenuation, ap[)[nepers/m], 
as 


2 


(5) 


rr 


a RD = 


aRpdr [nepers], 
o 


Several additional assumptions will be incorporated in the following 
development. These include the assumption that both the specific at- 
tenuation, a, and the reflectivity, Z, may be represented at any fre- 
quency by a power law dependence on the rain rate, R [rm/hr], i.e., 

b 

(6) a RD = a RD R 

RM RM 

and 

b z 

(7) Z = a z R 



where the subscripts RD and RM indicate that the subscripted quantity 
is to be evaluated at either the radar or radiometer frequency, re- 
spectively. Calculations based on Mie theory for various drop size 
distributions [4,5] indicate that this is a reasonable assumption 
even for frequencies where the drop sizes are resonant; this is dem- 
onstrated by Figures 1 and 2 which show regression power low fits to 
Mie calculations at frequencies from 1 to 500 GHz for Marshall -Palmer 
drop size distributions at 0°C. The a and b constants in Equations (6) 
and (7) are, of course, frequency dependent and have been tabulated for 
extensive frequency ranges in Reference 4 for the case of specific 
attenuation and in Reference 5 for both specific attenuation and 
reflectivity. 

It will also be assumed that attenuation and radiometric temperature, 
T s [°K], are related by 



where T m [°K] is the mean absorption temperature of the medium and the 
subscript » indicates that the attenuation is associated with a path 
extending to a resolution cell located at infinity [6]. This theoretical 
relationship is strictly valid for media where scattering is negligible; 
however, it has been demonstrated empirically that it is quite adequate 
for the rainfall case at frequencies at least up to 15 GHz [7,8]. An 
example of this agreement is shown in Figure 3 for the case of simul- 
taneous 15 GHz radiometric temperature measurements and attenuation 
measurements using the ATS-5 synchronous satellite beacon. 

Now, proceeding with the solution of the radar range equation, the 
rain rate, R, is eliminated from Equations (6) and (7) yielding 
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Figure 1. Specific attenuation, a , vs. rain rate 
at several frequencies. 
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Figure 2. Equivalent reflectivity, Z eq , vs. rain rate 
at several frequencies. 
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Figure 3. A comparison between measured attenuation and measured 
radiometric temperature at 15 GHz [8]. 
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(9) 


RD " a RD 


ft) 


RO 


Combining Equations (1), (5), and (9), the radar range equation Is 
expressed In terms of equivalent reflectivity: 


( 10 ) 


r P« « QZe 


. 2 ri«D f r Z b "D/»z dr 1 
.Iff. l», bRD/bz J ° J 


Making the substitutions 


. ?R0 

9 ' 2 az b RD/ b z 


- b RD/ b z 


(ID 

and 

(12) y * Z 
for simplicity, yields 

-b 2 /b RD 1 dr. 


(13) 


r 2 P r - Q y 


Now, the logarithm is taken of Equation (13); the result is then dif- 
ferentiated with respect to r and rearranged: 


(14) 


d y. . bop d[tn(r^P r )3 ^RO „ 
Hr + 6 r- • dr y b. 9 ' 


The solution of Equation (14) is 


b RD 

S7 


(15) 


c-^. g r <r*P r f*dr 


2n x b RD /b z 


(r%) 


where c is an unknown Integration constant. Substituting for g and y, 
the solution, „ 

r z P 

(16) Z • - r - 

IV.? a RD f r , 2„ , b RD/l>2 J b z /b R o 

[e 2 T7 -rr^L L (r p r) dp ] 


1 Sl b RD/ b z 0 


Is obtained. 
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The Integration constant, C, may be evaluated by considering the 
simple hypothetical reflectivity distribution shown in Figure 4. 



Figure 4. Reflectivity distribution used to evaluate the 
unknown constant of integration, c. 


Since the attenuation over the Interval from r=0 to r 0 Is zero, 
Equation (1) becomes 

(17) P r * % z o 
r o 

when evaluated for r=r 0 . Similarly, P r =0 for 0<r<r Q and, therefore, 
Equation (16) becomes 


(18) 


r 2p 

r o*r 


- b z/t>RD 


when evaluated for r=r 0 . Consequently, 
(19) c = Q b RD /b z . 


Substituting Equation (19) Into Equation (16), the following 
compact form of the solution of the range equation is obtained 


(20) 


z - r 


Q(l-I r ) bl/b| ® 
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where 


(21 ) 


2a RD t) RD 

M> I Q) l,RD/bz 



(r 2 P r ) bRB/bz 


dr . 


III. THE RADAR SATUR *ION FACTOR 

The quantity I r has a number of Interesting properties. For example, 
comparing Equations (1) and (20), it Is noted that 


( 22 ) 


A rd - m --i 


(1-Ir) 


b 2 /2b RD 


Thus, the quantity I r Is intimately related with the path attenuation. 
Furthermore, equating Equations (8) and (2z) for a resolution cell 
located at Infinity, It Is found that 

(23) 1 - li. (l-i ) bz/2bRD 

'm * 

where, again, the subscript » indicates that the integration Is to be 
performed over the entire range from rtO to *. Her? it is seen that 
I. is closely related to the ratio of the radiometric sky temperature, 
T s ,to the mean absorption temperature, T m , Since the ratio T $ /T m is 
a measure of the degree of saturation associated with the radiometric 
temperature measurement, it Is reasonable to associate that inter- 
pretation with the radar case as well. Thus, the quantity I r will be 
referred to as the radar saturation factor. 

The radar saturation factor, Ir, is unitless, range*- from zero to 
one, and Increases monotonically with range. Examining Equation [22) 
it will be noted that, as I r approaches unity, the path attenuation 
becomes extremely large. 

Examining Equation (21), it is noted that the calculation of Ir 
Involves a very simple range correction and integration of the received 
backscattered radar power. Thus, the real time calculation of I r 
provides a simple means of determining, in real time, whether a radar 
has become "blind" due to range attenuation and, furthermore, at 
what range the "blindness" occurs. 
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Now consider an extended region of uniform rain rate and, hence, 
uniform reflectivity as depicted in Figure 5. Examining Equation (20) 
It is clear tnat the quantity 


(24) 



I 

I 

l 

I 


is a constant independent of range since both Q and Z are independent 
of range. Furthermore, since the attenuation will increase. without 
bound as r increases, the range corrected received power, r^Pr. must 
approach zero (see Equation (1)). Therefore, the radar saturation 
factor must approach unity as the range increases. 



R 


Similar reasoning leads to the same conclusion even for the case 
of an inhomogeneous rain rate distribution if either the rain rate 
or extent is sufficiently large, i.e., 

(25) I r * i 

for any case where the path attenuation becomes sufficiently large. 

Substituting Equation (21) into (25) and rearranging yields 

(26) L. (yP r ) bRD/bz dr] VbRD * Q 


for sufficiently high or extended rain rates. The left side of 
Equation (26) has an absolute bound which is simply the calibration 
constant of the radar. Thus, if the left side of Equation (26) is 
evaluated and the maximum value found for experimental data obtained 
using an attenuating wavelength radar viewing a sufficiently intense 
storm, this value establishes an absolute radar calibration constant 
independent of any other measurements. 


IV. RADIOMETRIC TEMPERATURE 


Since the reflectivity, Z, can be found as c junction of range. 
Equation (20), the rain rate, R, is also known, tquation (7): 


rM /bz 

* a zQ J (l-I r ) 1/b RD 


Similarly, using Equation (6) the specific attenuation can be found 
at any other frequency. In particular, let the specific attenuation 
be evaluated for some radiometer frequency other than the radar fre- 
quency; then. 


(28) 


a RM = a RM 



1 

b RM/ b RD 

-I r ) 


and the path attenuation at the radiometer frequency is 


(29) 


. a RM 

RM » " (a z Q) b RM /b z 


O b RM/bz 
(r P r ) 


(i-i r ) bRM/DRD 


dr. 
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Combining Equations (8) and (29) yields 


(30) 


T s a T^l-exp 


(a z Q) 


b RM /b z 


(r 2 P r ) 

(1-I r ) 


b RM /b z 

b RM^ b R0 


v 

so that the radiometric sky temperature can be obtained 
directly from the radar data. 



at any frequency 


V. RADAR CALIBRATION 


If Equation (30) is rearranged, the radar calibration constant, 

Q, can be expressed in terms of the backscattered radar power and the 
radiometric temperature as: 


(31) 



a RM 


r <r2p r ) bRH/bz 

0 (l-I r ) b RM /b RD 



dr 


b 2 /bRM 


The right hand side of this expression is implicitly dependent upon Q 
because of the presence of Ip. Nevertheless, if a reasonable estimate 
of Q can be made, this expression can be used to solve for Q in a 
stable, iterative manner. The estimate of Q should be used to cal- 
culate the radar saturation factor, lr, and then the improved estimate 
of Q is calculated using Equation (31). Thus, a bore-sight radiometer 
operating at a frequency either the same as or different than the radar 
frequency can be used to calibrate the radar on an instantaneous basis. 

Obviously, if the radar is operating at a non-attenuating wavelength 
or If the rainfall observed is very light, the radar saturation may 
be very small or even negligible; in this case, neglecting I r in 
Equation (31) eliminates the need for an iterative solution of this 
equation. 

Conversely, if very intense rain is observed at an attenuating 
wavelength, the radar saturation factor, I r , will approach unity for 
large r; in this case the Integrand in Equation (31) will become in- 
determinate and, thus, may generate significant errors. However, 
this Is precisely the case where the bound given by Equation (26) is 
applicable; thus, the radiometric temperature measurement and the 
use of Equation (31) are not required. 
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VI. SUMMARY 


Using the Hitschfeld-Borden approach, the radar reflectivity, Z, at 
a given range has been obtained in terms of the received radar power, 
the radar calibration constant, and the parameters of the power law 
relations between specific attenuation, reflectivity, and rain rate. 

This solution of the radar range equation contains a factor, I r , re- 
ferred to as the radar saturation factor, which is dependent upon the 
integral of the range corrected received radar power. This radar 
saturation factor, in turn, is shown to exhibit many interesting char- 
acteristics; for example, this factor is very simply related to the path 
attenuation and the radiometric temperature. Furthermore, this factor 
is shown to saturate, i.e., approach a known bound, in a manner analo- 
gous to that in which radiometric temperature saturates. Finally, these 
properties are used to develop two alternative techniques useful for the 
determination of the radar calibration constant. The first technique, 
applicable to radars operating at attenuating wavelengths, permits a 
direct, absolute calibration of the radar in terms of the measured radar 
data itself. The second technique, employing independent radiometric 
temperature measurements, permits continuous radar calibration during 
data acquisition. 
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